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1.  INTRODUCTION 

The  days  of  laborious  hand  extraction 
of  particle  size  distribution  Information,  partl- 
by  particle,  arc  largely  over.  Now  that  autceiatlc 
equipment  "snows"  us  with  numbers  we  must  learn  - 
or  teach  the  ccmputer  - how  to  extract  swanlngful 
values  from  this  avalanche  of  data.  This  paper 
rcTports  on  efforts  in  this  direction. 

With  the  advent  of  such  copious 
particle  data  it  is  becoming  more  realistic  to 
plan  for  and  Obtain  Independent  measurements  of 
related  parasieters  both  from  aircraft  and  from 
remote  sensors  on  the  ground.  Some  comparisons 
of  this  nature  are  given. 

2.  DATA  SOURCE 

For  several  years  the  AFGL  has  flown 
a full  set  of  Particle  Measuring  Systems,  Inc. 

(PMS)  sensors  aboard  both  a C-130  (AF  Operated) 
and  a Learjet  (contractor  operated) . Many 
flights  have  been  made  through  U.S.  east  coast 
deep  cyclonic  systems  as  well  as  Into  tropical 
cloud  systems.  A variety  of  hydrometeor  types 
has  been  encountered.  Many  of  these  aircraft 
data  samples  were  accompanied  by  measurements  of 
the  radar  reflectivity  factor  just  ahead  or 
behind  the  aircraft  made  by  a ground  based  radar, 
as  well  as  measurements  of  the  water  content  and 
particles  by  other  airborne  instrumentation. 


3.  THE  PMS  OPTICAL  ARRAY  PROBES 

The  PMS  optical  array  instruments, 
and  their  calibration  and  to  some  extent  data 
analysis  have  bean  described  and  discussad  In  a 
number  of  publications,  reports  and  preprints, 
such  as  Knollenbarg  (1970) (1972) (1975)  (1976) , 
Haymsfiald  and  Knollenbarg  (1972) , Raymaflald 
(1976) , Cunningham  (1976) , Gayat  (1976) . Aspects 
of  tha  design  which  relate  to  the  analysis 
procedures  described  in  this  paper  follow. 

A vortical  laser  beam  llluminatas 
through  an  optical  system  a line  of  photodiodes, 
24  in  tha  case  of  AFOL's  ID  probes  and  32  for  tha 
20  probes.  Particles  on  passing  through  this 


beam  cast  a shadow  on  these  diodes.  If  the 
shadow  is  complete  enough  to  reduce  the  output 
level  of  the  diode  circuit  by  1/2  then  this 
diode  Is  considered  occluded.  In  one  set  of 
Instruments,  the  ID's,  the  maximum  number  of 
diodes  occluded  during  the  passage  of  a particle 
Is  recorded  as  particle  size.  In  the  other  set, 
the  2D's,  a rapid  scanning  system  is  used  to 
allow  the  recording  of  the  diodes  occluded  per 
2Sum  (cloud  probe)  or  200)im  (preclp.  probe) 
forward  motion  of  the  aircraft.  A plctograph  of 
the  particle  can  be  produced  after  suitable  data 
processing  of  this  latter  type  of  information. 

The  end  diodes  in  ID  devices  are  used  to  indicate 
oversize  particles  which  permits, through  suit- 
able design,  the  electronic  hardware  to  reject 
particles  which  are  only  partly  saiRpled.  This 
same  problem  Is  handled,  in  the  case  of  the  2D 
devices,  by  a subsequent  software  routine. 
Unfortunately,  with  this  edge  reject  requlresient 
sample  volume  decreases  with  Increasing  particle 
size.  Schemes  can  be  devised,  which  so  far 
require  visual  examination  of  the  particle, 
that  could  estimate  particle  size  from  the 
plctograph  of  a portion  of  the  particle  and 
thus  retain  the  full  sasqile  volume  for  the  larger 
fewer  particles.  It  has  not  been  found  practical, 
particularly  for  Irregular  particles,  to  pursue 
this  approach. 

The  three  length  dimensions  that  make 
up  the  sample  volume  are  detenilned  by;  one,  the 
length  of  the  diode  array  divided  by  the  magnifi- 
cation euid  adjusted  for  the  and  reject  feature; 
two,  the  distance  of  eMPOsurs  of  the  laser  beam 
to  the  free  flowing  air  and  particle  stream,  ad- 
justed for  a depth  of  field  factor  for  the  small 
particles  and;  three  the  length  correaiwnding  to 
the  chosen  sample  time  multiplied  by  true  air- 
speed. Length  one  and  three  are  determined  with 
ease  as  long  as  the  isagnlfl cation  of  tha  optics 
Is  checked  frequently  and  good  true  airspeed 
values  are  available.  The  second  length  has  been 
tha  subject  of  soma  debate  because  of  the 
requirement  in  the  smaller  channel  sizes  for 
depth  of  field  Infozziatlon.  Irregular  particles 
with  thin  arms  or  particles  with  large  aspect 
rations  may  present  a problem.  In  eur  present 
analysis  we  have  used  tha  depth  of  field  values 
given  by  the  manufacturer,  PMS. 
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4.  ANALYSIS  OF  ID  DATA 

4.1  Adjufltnent  for  Qu«ntl»«tlon  and  for 

Non  Spherical  Sh«p«« 

Several  analytical  steps  axe  required 
before  the  counts  per  channel  output  of  the  PMS 
circuitry  can  be  processed  Into  neaningful 
particle  distribution  Information.  First  the 
basic  magnification  of  the  optical  Imaging 
system  must  be  known.  This  Is  frequently  check- 
ed, and  adjusted  If  needed,  by  passing  a length 
of  wire  of  known  size  through  the  beam.  Because 
the  particle  size  Information  Is  derived  from 
discrete  steps  and  because  there  is  some  lost 
space  between  photodiodes  there  Is  a quantization 
correction  to  be  applied  even  In  the  case  of 
spherical  particles  (rain) . in  the  case  of  non 
spherical  shapes.  Irregular  or  regular,  further 
corrections  are  needed.  A study  of  the  correc- 
tions required  to  allow  conversion  from  irecorded 
channel  (diode)  number  to  maximum  particle  size 
was  supported  by  AFGL  and  made  t>y  Knollenberg 
(1975) . The  conversion  equations  for  the 
particle  types  used  In  the  AFGL  storm  flights 
have  been  derived  In  part  from  this  study  and 
Incorporated  Into  our  software.  Types  "Plata 
family"  and  Aggregates  of  Plates  and  dendrites 
have  recently  been  added  to  the  relations  listed 
for  thesa  being  estimates  using  the  standard 
more  pristine  types  as  a guide.  The  equations 
being  used  by  AFGL-LY  are  listed  In  Table  1. 

Table  1 

Adjustment  of  PMB  Channel  Number 


Adjusted  Ch  t • 

* Ind. 

Ch  # ♦ 

B 

Breakpoint 

Particle  ^pe 

M 

B 

(Ind. 

Ch  •) 

Rain 

.99 

.18 

- 

Net  Snow 

.99 

.18 

LE 

2 

N H 

1.15 

.18 

GT 

2 

Large  Snow,  Small  Snow 

1.15 

.18 

- 

Bullet-Hosattes 

1.02 

.32 

- 

Columns  (4il) 

1.3 

.76 

- 

Needles  (7.5:1) 

.20 

3.04 

LE 

1 

N « 

1.28 

1.08 

GT 

1 

Plata  Family 

.94 

.55 

LE 

3 

H N 

i.oe 

.02 

or 

3 

Agg.  Plates  and 

Dendrites 

.94 

.55 

LB 

3 

m m n n 

1.20 

-.44 

GT 

3 

4.2  Conversion  of  Crystal  Size  to 

Eqt>lvalent  Malted  Diameter 

Use  of  the  above  relations  permits 
analysis  of  the  "cloud"  and  "preclp"  PMS  probe 
raw  counts  to  detazmlna  particle  spectra  In 
for  Instance,  the  form  of  maximum  physical  size 
vs  nimiber  of  particles  per  cu  meter.  In  order 
to  derive  normalised  particle  spectrum,  mass 
spectrum  or  total  Ice  or  liquid  water  content 


(INC,  LHC)  or  radar  reflectivity  (Z)  an  addition- 
al conversion  Is  required,  that  of  maximum  size 
to  mass  or  more  conveniently  to  equivalent  malted 
diameter  (D) . 

There  Is  now  a considerable  literature 
reporting  on  such  relations.  Unfortunately,  the 
data  set  Is  somewhat  biased  toward  single  prlsteen 
crystal  measurements  and  rimned  particles.  It  Is 
also  geographlcetlly  biased  toward  mountain  snow 
systems  and  against  deep  cyclonic  storm  cases.  A 
listing  of  the  several  sources  and  Illustrations 
showing  some  of  the  relations  found  Isetween  maxi- 
mum length  and  particle  mass  has  recently  been 
given  by  Heymsfleld  (1976). 

From  the  above  noted  sources,  and  a 
few  others,  a set  of  power  functions  has  been 
derived.  Those  found  most  useful  In  our  analysis 
of  east  coast  deep  storm  clouds  are  given  In  Table 
2.  niree  rather  different  particle  type  L to  D 
relations  eure  also  shown  in  Fig.  1 along  with 
marks  denoting  the  original  channel  center  size 
and  adjusted  size  for  channels  1,  5 and  IS. 

Table  2 

Equations  used  to  Convert  Particle  Size  (L) 


to  Halted  Diameter 

(D),  D 

w AL^  Units,  am. 

Particle  Type 

A 

B 

Breakpoint  (L) 

Rain 

1 

1 

- 

Net  Snow 

1 

1 

IZ  1 

H *1 

1 

.65 

GT  1 

Large  Snow 

.4 

.78 

IZ  1 

(1  tf 

.4 

.88 

GT  1 

Small  Snow 

.4 

.78 

IZ.5 

n n 

.37 

.67 

GT.5 

Bullet-Rosettes 

.26 

.67 

LB. 2 

M n 

.44 

1.0 

GT.2 

Columns  (4:1) 

.44 

1.0 

- 

Needles  (7.5:1) 

.26 

.67 

- 

Plate  Family 

.34 

.78 

LE  1 

N N 

Agg.  Plates  and 

.34 

.68 

GT  1 

Dendrltea 

.34 

.78 

LE  1 

ti  n tt  ft 

.34 

.68 

GT  1 

Two 

rather 

broadly  defined  snow 

classes  not  found  explicitly  In  standard 
classifications,  l.e. , Magono  and  Lae  (1966) 
are  large  snow  (LS)  and  small  snow  (SS) . It  was 
found  In  practice,  particularly  whan  probing 
deep  cyclonic  systems  that  these  types  bast 
dascrUfte  the  snow  and  snow  aggregates  that  could 
be  Identified  from  a rapidly  isovlng  Instrumented 
aircraft.  This  slnpllficatlcn  Is  necessary  only 
In  part  because  of  the  nature  of  the  observation 
platform..  When  dendrites,  plates  or  needles 
were  encountered  they  could  In  fact  be  identified 
from  the  aircraft  on  the  so  called  "snow  stick" 
or  on  the  2D  real  time  CRT  display.  It  had  also 
been  found  In  our  work  with  INC  and  radar  Z 
correlations to  bo  useful  to  lusg)  several  of  the 
available  relations  Into  these  two  categories. 
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ANALYSIS  OF  2D  DATA 


riguA-t  ).  Ploti  oi  convtJuion  Aelatiotu  b^XMeen 
poAticZe.  ilzz  ID  a>id  t^eJUUd  VhmvUA  (P)  OiUgAtml 
chatm&l  iAze.  muttLptie.d  by  diode  uiixLth  ii 
indicated  on  L axli.  CP,  cloud  pAobe,  PPj  pAecip 
pAobe;  PF,  plate  iamily;  N,  needlei. 


The  class  LE  Is  loosely  defined  as 
aggregates  of  nany  plates  and  undefined  crystals. 
Large  snow  clusters  (cm.  In  slse)  often  found  In 
and  just  above  the  melting  layer  (In  the  radar 
'bright  band')  probably  should  not  be  classed  as 
LS,  but  deserve  a class  by  themselves,  this  has 
not  been  defined  as  yet.  One  reason  for  suggest- 
ing a separate  class  of  large  snow  clusters  Is 
that  It  would  probably  have  a Z to  M relation 
significantly  different  than  that  used  for  LS  or 
single  crystal  types.  The  class  SS  Is  also 
loosely  described  as  undefined  single  crystals 
or  aggregates  made  up  of  only  a few  single 
crystals. 


4.3  ID  Output  Format 

nie  one  second  count  In  each  of  the 
15  channels  of  each  probe  Is  suiMd  for  a 
preselected  time  period.  Four  seconds  Is 
usually  chosen  for  this  time  period  as  this  Is 
roughly  the  time  It  taXea  the  aircraft  to  fly 
through  one  sample  volume  of  the  ground  radar. 
Based  on  calculations  of  Joss  and  Naldvogel 
(1969)  the  sasg>le  volume  swept  out  In  four  sec. 
by  the  PMS  optical  array  probes  Is  not  a suffi- 
cient sample  size  for  the  larger  particles  of 
each  sensors  range.  In  order  to  obtain  an 
improved  distribution  function  we  have  usually 
combined  the  four  sec.  samples  Into  5 min.  "run" 
samples  normalizing  the  four  sec.  data  In  a 
manner  suggested  by  Sakhon  and  Srlvastava  (1970) . 
The  results  of  such  a treatment  have  been 
Illustrated  by  Cunningham  (1976) . It  should  be 
noted  that  the  coa^utad  water  content  and  even 
radar  reflectivity  values  froei  the  four  sec. 
sample  appear,  hcwever,  to  give  meaningful  and 
not  particularly  statistically  noisy  values. 

Sea  Fig.  1,  Cunningham  (1976). 


As  discussed  above  the  2D  data  format 
contains  Information  on  the  shadow  of  the  particle 
In  two  dimensions,  one  along  the  effective  plwe 
of  the  diodes  the  other  In  the  direction  of 
flight.  A number  of  types  of  output  can  be 
derived  from  this  Information,  the  type  usually 
presented  Is  the  direct  plctograph  of  the 
particle  shadow.  He  have  chosen  to  use  this 
presentation  and  several  others.  Before  de- 
riving the  more  sophisticated  parameters  eui 
edit  program  has  been  found  to  be  a necessity. 


5.1  2D  Edit  Program 

The  plctograph  output  is  very  useful 
%d.thout  any  editing  of  the  data  as  one's  'eye' 

Is  very  generous  In  its  ability  to  automatically 
throw  out  extraneous  noise,  malfisictions,  or 
usually  unwanted  effects  from,  for  Instemce, 
partly  transparant  crystals.  Programs  have 
been  written  to  Indicate  the  occurrence  of  mal- 
functioning diodes,  to  eliminate  particles  when 
splashing  or  breaJcup  occur,  and  to  fill  central 
holes  In  particle  images.  This  becomes  a problem 
with  rain  data  as  well  as  with  transparent  ice 
particles  as  rain  drops  transmit  light  which 
often  results  In  an  apparent  hole  In  the 
particle  plctograph. 

Several  other  'fixes'  are  necessary. 
Tlaia  Is  required  for  triggering  the  scemnlng 
circuits  as  the  particle  begins  to  shadow  the 
diode  array.  Thus  part  of  the  particle  is  not 
Included  In  the  final  plctograph.  A correction 
has  been  added  In  the  software  that  cooqpensates 
for  this  missing  part.  Area  is  added  as  a 
fraction  of  the  number  of  diodes  shadomd  on 
the  first  real  scan,  alternately  adding  one  scan 
or  two  scans  to  each  successive  frame. 


At  times  drops  hit  some  of  the 
sensor  structure  and  then  enter  the  Illuminated 
sample  region.  These  drops  appear  as  long 
strea]cs.  They  are  partly  eliminated  in  the 
analysis  by  rejecting  all  particles  larger  than 
32  diodes  In  the  'time*  direction. 

Finally  the  equivalent  of  the  end 
diode  reject  procedure,  built  Into  the  ID 
hardware.  Is  in  this  case  Included  In  the  2D 
software. 


5.2 


Size  Determined  In  Three  Wavs 


SMtta 


In  our  present  procedures  three 
particle  size  parameters  are  computed;  one,  the 
equivalent  of  the  ID  alze;  two,  an  approximation 
of  the  sHudmum  length  of  the  particle  shadow  on 
the  horizontal  plane  and;  three,  the  area  of  the 
shadow  In  the  horizontal  plane.  The  maximum 
length  Is  approximated  by  the  following  ej^reasion 


□ D 


(X  +v'X®  + y‘')/2  it  X>Y  ; 


(Y+yx*+Y*)/2  if  Y>X 


where  X ia  the  maximum  number  of  diodes  shadowed 
in  the  horizontal  direction  and  Y is  the  meuci- 
mum  number  of  scans  shadowed  in  the  flight 
direction.  Analytical  experiments  with  more 
sophisticated  derivations  of  a maximum  length  of 
irregular  shapes  show  that  this  simple  formula 
ia  quite  adequate.  If  the  particle  type  is 
noted  as  rain, just  the  average  of  X and  Y is  used. 


The  area  shadowed  by  the  particle  is 
conveniently  obtained  by  just  scanning  the  diodes 
occluded  on  each  scan,  after  the  editing  programs 
mentioned  above  have  been  applied,  then  determining 
the  total  sum. 


12  3 4 

E.C.R. 


Derived  Shape  Parameters 


Two  shape  parameters  have  been  under 
study.  One  is  a measure  of  edge  cosplexlty  of 
the  particle , 1 . e . , the  ratio  of  the  average 
projected  size  of  the  pictograph  vs  its  perimeter 
has  been  named  average  projection  ratio.  The 
other  la  called  the  equivalent  circle  ratio  (ECR) 
and  is  the  ratio  of  the  perimeter  of  the  picto- 
graph to  its  area  divided  by  the  ratio  of  the 
perimeter  of  a circle  to  its  area,  the  areas 
being  equal.  These  ratios  are  coaputed  for  each 
sample  run  (about  5 minutes)  and  tabulated 
against  the  meuclmum  size  parameter.  These  tables 
are  being  studied  to  determine  the  feasibility 
of  using  such  shape  functions  to  determine  snow 
type  automatically.  Snow  particle  shapes, 
modeled  after  the  patterns  used  by  Knollenberg 
(1975)  were  made  for  several  sizes,  relative  to 
the  32  diode  2D  dimension,  and  the  above  ratios 
determined  from  a 32x32  matrix  of  squares  for 
various  equivalent  maximum  sizes  and  orientations 
of  the  particle.  The  more  significant  ratio 
appears  to  be  the  ECR.  Curves  of  its  value  vs 
maximum  size  are  shown  in  Fig.  2.  An  example 
of  2D  data  superimposed  on  these  dashed  lines 
indicates  that  in  this  case  the  preponderance 
of  snow  type  was  LS  and  6S. 


5.4  JUrea  to  Malted  Diameter 

Particle  area,  or  at  least  the  pro- 
jected horizontal  area,  can  as  has  bean  mention- 
ed, be  calculated  using  the  2D  output.  Calcula- 
tions, independent  of  our  previous  L to  D 
conversion,  can  then  Im  made  to  determine  particle 
mass  or  malted  diameter.  Relations  iMtween 
particle  area  and  mass  are  required  as  very  few 
were  available  Soma  were  obtained  by  rafozmating 
wor)(  by  Auer  (1971) . New  data  were  obtained  from 
Knollenberg  (1976) . From  these  two  sources  a 
selection  of  particle  types  was  made  that 
showed  distinct  diffaranoas.  Soma  additional 
functions  were  added  for  types  Which  exist  over  a 
limited  else  span,  i.a.,  wet  snow,  through 
reference  to  ttelr  L to  D relations.  The  result- 
ing relations  for  these  selected  types  are  given 
in  Table  3. 


FiauA.e  2.  Pio-t  zquivatetit  cJjuJLt  >uvUx>  (ECR) 
n'OMjnm  i^ize.  ioA.  nonujiat  aujitxiJL 

ihapzi:  C,  coturna;  R,  AOieXtu;  LS,  SS,  LoAnz 
tuid  imatt  4Kow;  SP,  izctoA  platz;  PP,  pintiti 

d.et'.dAitz. 


Table  3 

Equations  used  to  Convert  Particle  Area  (A) 
to  Melted  T^lametar  (D) , D - aA^  Units;  mm,  sn^ 


Particle  Type 

a 

b 

Brea)cpolnt 

Rain 

1.13 

.5 

- 

Net  Snow 

1.13 

5 

LE.78 

tf  tt 

1.08 

.32 

Between 

H n 

.52 

.40 

or  62 

Large  and  Ssiall  Snow 

.52 

.40 

GT.OOl 

Bullat-Rosettes 

.50 

.39 

LE.21 

M m 

.61 

.51 

0r.21 

Columns 

.42 

.35 

U.16 

H 

.52 

.46 

GT.  16 

Needles 

.42 

.35 

- 

Plata  Fasd.ly 

.46 

.32 

LE.45 

M m 

Agg.  Plates  and 

.56 

.58 

01.45 

Dendrites 

.39 

.31 

LB.  25 

N It  n n 

.44 

.40 

01.25 

Dendrite  Family 

.39 

.31 

LF  1.07 

t<  tl 

.38 

.48 

CT  1.07 

6.  COMPARISON  OF  RESULTS 


It  is  always  of  interest  as  well  as  a 
certain  measure  of  progress,  or  lade  thereof,  to 
compare  results  between  different  ways  of  sensing 
a quantity  or  of  analyzing  data.  A few  axasples 
of  such  comparisons  are  given  below  from  a vary 
limited  set  of  available  data. 
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1 
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1009 


TIME  (?) 


Tlc.uJie.  j.  ComfxvuJton  te-iit'eei:  Z"'C  c^e.'Uv&d 
P!'S  poAtictt  data  and  that  ifiom  thz  H'lEf! 
■uiit/iunznt. 


pariaon  batwean  tha  ID,  2D  Input 


Coopaxlaon  of  results  batwean  the 
sevaral  ways  of  anslysls  outllnad  abova  can  ba 
made,  hn  intazcospaxlaon  batwaan  the  neaaured 
total  particle  count  cooiputed  INC  (H)  and  radar 
raflactlvlty  (Z)  is  givan  In  table  4,  The 
results  fm  two  ccctplataly  separate  Instrunents 
are  shorn  hera  tha,  ID  and  211  as  well  as  three 
diffarant  ways  of  coi^nitlng  M and  Z from  the  2D 
device. 


near  Nallops  Island  VA  illustrates  both  fair 
agreement  imd  the  necessity  to  consider  the  con- 
tributions of  particles  larger  than  the  probe 
can  properly  sai^la.  Table  S compares  the  M and 
Z values  computed  from  the  ID  device  as  well  as 
values  obtained  for  an  axtri^latad  distribution 
(m  , Zco) , distribution  Is  extrapolated 

using  as  a guide  the  shape  of  the  nondlmenslonal 
distribution  function.  The  remotely  measured 
values  of  radar  reflectivity  (Z^^)  were  obtained 
from  a carefully  calibrated  ground  radar.  Note 
the  Z values  lie  as  they  should  between  Z and 
Zoo  on  TOst  of  the  passes. 


6.2  Comparison  of  ID  M and  Z Results  with 

Nmmote  Radar  Msasureisents 

An  axaspla  of  good  agreement  between 
PMS  derived  Z values  and  ground  radar  values  are 
given  by  Cunningham  (1976) . The  agreement  was 
better  In  the  matdilng  of  the  details  of  Z struc- 
ture than  perhaps  In  absolute  value,  however,  the 
difference  converted, for  Instance,  to  L to  D 
error  would  be  small.  Cosparlson  between  derived 
Z and  radar  measured  Z made  on  the  20  March  77 


Table  4 

comparison  of  Four  Derivations  of 
Paramatars  of  a Particle  Distribution 

(Data  Taken  from  20  March  77  Flight  of  C-llO 
10;ia;22  - 10:18:57) 


T 3 
#/n 

"3 

q/m 

dBZ 

ID 

1800 

.97 

36 

2D,  HFP*^’ 

2220 

.93 

35 

2D,  Max.  Length 

2417 

1.25 

37 

2D,  Area 

7.288 

.98 

35 

Cable  i 

C-130  Data  from  20  March  77  Dallops  Is  Va. 
Snow  Type  Aggregates  of  Plates  and  Dendrites 


Start 

Time 

Z 

Alt 

Km 

Temp 

oc 

M 

m 

(1) 

dBZ(l) 

dBZ^ 

1006 

6 

-17 

.22  - 

.25 

17 

- 

24 

23 

0951 

4.9 

-13 

.37  - 

.39 

22 

- 

24 

24 

0947 

4.9 

-13 

.23  - 

.23 

21 

- 

27 

24 

1015 

4.9 

-13 

.27  - 

.39 

22 

- 

31 

33 

0937 

4.1 

-10 

.19  - 

.23 

21 

- 

25 

25 

0933 

3.9 

- 9 

.23  - 

.27 

22 

- 

25 

26 

1028 

3.9 

- 9 

.24  - 

.30 

22 

- 

27 

27 

0917 

3.0 

- 5 

.08  - 

.26 

19 

- 

35 

29 

(1)  Data  oroceased  In  a slnllar  way  as  ID  Cate. 


(1)  First  number  value  calculated  from  PMS 
data,  second  number  calculated  after 
extrapolating  distribution  to  Infinite 
particle  site. 

(2)  value  measured  bv  remote  ground  based 
radar,  SPANDAA. 
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